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In the presence of time- reversal symmetry, quantum interference gives strong corrections to the 
electric conductivity of disordered systems. The self- interference of an electron wavefunction trav- 
eling time-reversed paths leads to effects such as weak localization and universal conductance fluc- 
tuations. Here, we investigate the effects of broken time-reversal symmetry in a driven artificial 
two-level system. Using a superconducting flux qubit, we implement scattering events as multi- 
ple Landau- Zener transitions by driving the qubit periodically back and forth through an avoided 
crossing. Interference between different qubit trajectories give rise to a speckle pattern in the qubit 
transition rate, similar to the interference patterns created when coherent light is scattered off a 
disordered potential. Since the scattering events are imposed by the driving protocol, we can control 
the time-reversal symmetry of the system by making the drive waveform symmetric or asymmetric 
in time. We find that the fluctuations of the transition rate exhibit a sharp peak when the drive 
is time-symmetric, similar to universal conductance fluctuations in electronic transport through 
mesoscopic systems. 



At low temperatures, the conductivity of disordered 
systems is strongly influenced by quantum interference 
effects such as weak localization (WL) and universal con- 
ductance fluctuations (UCF) p~H3]. Weak localization is 
due to constructive self-interference of an electron wave- 
function traversing time-reversed paths, leading to an en- 
hanced probability of back-scattering and therefore a re- 
duction of the conductance. UCF describe the strong 
fluctuations in conductance that occur as a function of 
any parameter that changes the scattering configuration 
[3H6] . The interference effects are highly sensitive to any- 
thing that breaks time-reversal symmetry, such as a mag- 
netic field applied perpendicularly to the motion of the 
charge carriers. Studies of weak localization and UCF 
thus provide a method for investigating effects related to 
phase coherence, which has been used in a wide variety 
of systems ranging from metals [7] and semiconductors 
[8] to quantum dots [9HTT] and graphene [12], and even 
for the scattering of light off disordered media [T3UT5] . 

The presence of weak localization and UCF in such 
widely varying systems shows the universality of the ef- 
fect, occurring independently of the sample size, dimen- 
sionality and the degree of disorder. In mesoscopic sys- 
tems, there are typically a large number of scatterers, 
giving millions of interfering paths that contribute to the 
electron transport. In this work, we investigate the effect 
of time-reversal symmetry on a driven artificial two-level 
system in the few-scatterer limit. The scattering events 
are implemented as Landau- Zener transitions by driv- 
ing a qubit multiple times through an avoided crossing. 
With both the number of scattering events and the time- 
reversal symmetry imposed by the driving protocol, we 
have control over the number of possible paths in the 
system. In a configuration with only four scatterers, we 
measure a sharp increase in the fluctuations of the qubit 
transition rate when the drive waveform is made sym- 



metric in time. 

We use a superconducting flux qubit [16] H7], which 
consists of a niobium loop interrupted by three Josephson 
junctions [18 , with a magnetic flux $ threading the loop. 
The clockwise and counterclockwise persistent currents 
±ip, corresponding to the qubit's diabatic states, are 
tunnel coupled with strength A. The two- level Hamilto- 
nian H = (— 1/2) (Aa x + sa z ) describes the qubit dynam- 
ics, where e = 2Ipf, f = 3> — <£o/2 is the flux detuning, 
$o = h/2e the superconducting flux quantum, and <r x 
and <j z are Pauli matrices. The qubit energy separation 
£?oi = V £ 2 + A 2 is therefore controlled by the flux <l> in 
the loop [Fig.Qa)]. 

We first discuss the concept of scattering events in a 
driven two-level system. Starting at negative flux detun- 
ing and with the qubit in its ground state, we apply a 
large-amplitude flux signal that drives the qubit through 
an avoided crossing and back again [Fig. [lja)]. At the 
first avoided crossing, the ground state |0) undergoes a 
Landau-Zener transition and splits into a coherent super- 
position of |0) and |1). The states evolve independently 
until the second time they reach the avoided crossing, 
where they interfere constructively or destructively de- 
pending on the relative phase cp acquired between the 
two transitions. 

The Landau-Zener transitions and the qubit evolution 
can be thought of as a phase-space analog of an optical 
Mach-Zehnder interferometer [18]. Figure [TJb) shows the 
energy evolution of the qubit during the drive, where 
the interference phase (p = (1/K) J^ 2 Eoi(t) dt is given by 
the shaded area between the two scattering events. The 
setup is similar to the mesoscopic system of Fig. [TJc) in 
the limit of very few scatterers. The qubit phase space 
of Fig. [TJb) contains only two scatterers, two possible 
trajectories and one interference phase. The problem can 
be solved analytically, with the resulting qubit transition 
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FIG. 1: Qubit trajectories during driven evolution. (a) 
Energy- level diagram of a flux qubit. Driving the qubit 
through the avoided crossing induces Landau-Zener transi- 
tions between the two states. The states recombine when the 
qubit is brought back through the avoided crossing, with the 
final outcome depending on the phase accumulated during the 
flux excursion, (b) Qubit transitions visualized as two scat- 
tering events in a Mach-Zehnder interferometer setup, (c) 
Illustration of a mesoscopic system with two paths related by 
time-reversal symmetry, (d-e) Qubit transitions for a drive 
waveform that brings the qubit through the avoided cross- 
ing multiple times. The blue and red traces mark two (out 
of eight) possible trajectories in the system, (d) When the 
waveform is symmetric in time, the trajectories acquire the 
same phases (since ipi = ^3), and will interfere constructively, 
(e) Without time-reversal symmetry, (pi ^ ips. 



rate showing oscillations as a function of the interference 
phase (f [T9] . 

We can increase the number of trajectories in our sys- 
tem by driving the qubit back and forth through the 
avoided crossing several times for each cycle of the driv- 
ing waveform. This allows us to increase the complexity 
of the system and move closer to the mesoscopic case 
while still having a controlled, non-chaotic phase space. 
Figure [TJd) shows the qubit being driven through the 
avoided crossing two additional times, giving a total of 
four scatters, three interference phases, and 2 3 = 8 possi- 
ble paths. Note that the drive waveform is symmetric in 
time. Qubit trajectories that are related by the temporal 
symmetry will pick up the same phase during the driven 
evolution, and they will therefore interfere constructively. 
Examples of time-reversed paths are plotted in red and 
blue in Fig.[jjd). We can easily break the time-reversal 
symmetry by making the drive waveform asymmetric in 
time, as shown in Fig.[TJe): the red and blue trajectories 
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FIG. 2: Qubit transition rate as a function of the time-reversal 
symmetry of the drive waveform, (a) Drive waveform for dif- 
ferent values of the drive asymmetry parameter a. (b) Zoom- 
in around the center region of (a); the waveform becomes 
asymmetric in time when a 7^ 0. (c) Measured transition 
rate versus static flux detuning fdc and a. (d) Qubit energy 
diagram, (e-f) Case when fdc > (positions in (c) marked 
by I and II) ; the qubit is driven through the avoided crossing 
twice, resulting in Landau-Zener oscillations as a function of 
the accumulated phase (p. (g) Case when fdc < (position 
III); here, the qubit is driven through the avoided crossing 
four times, giving rise to an intricate interference pattern. 



will acquire different phases since cpi and cps are no longer 
equal. The qubit thus provides a well- controlled test sys- 
tem for investigating the effects of broken time-reversal 
symmetry in the few-scatterer limit. 

We implement the drive protocol using a biharmonic 
signal [20], 

f(t) = fdc + M cos(ut + a) - A 2 cos(2ut). (1) 

We fix the frequency uj/(2tt) = 125 MHz, the amplitude 
Ai = 3m$ and the amplitude ratio A 2 /Ai = 0.55, 
while the parameter a controls the waveform's asym- 
metry. The waveform is plotted in Fig. |5|a) for a few 
different values of a, with Fig. [2^b) showing a magnifi- 
cation around the time t = 0. Note that the function is 
symmetric in time for a = and becomes increasingly 
asymmetric as a is increased. 

The waveform's period 2tt/uj = 8 ns is comparable to 
the qubit's dephasing time w 10 ns. This puts the dy- 
namics in the quasiclassical regime [19 , where coherence 
is preserved within one period of the drive, but where 
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multi-photon processes due to coherence over many pe- 
riods are not resolved [T8] . However, since the energy- 
relaxation time T\ « 20 [is is much longer, consecutive 
periods of the drive signal will lead to a build-up of 
excited-state population. The process of reaching the 
equilibrium population follows an exponential time de- 
pendence, which we characterize by a transition rate T. 
In the relevant regime (T ^> 1/Xi), we have T = 2W, 
where W is the probability (per unit time) of inducing 
a qubit transition from |0) to |1) after one period of the 
drive waveform [T9] . 

The measurement procedure consists of three steps: 
first, the qubit is cooled to its ground state using a 3- 
/is cooling pulse [21]; next, we apply the drive signal to 
induce qubit transitions; finally, the qubit state is read 
out using a superconducting quantum interference de- 
vice. By changing the length of the drive pulse and fit- 
ting the measured qubit population versus pulse duration 
to an exponential decay, we extract the transition rate T. 
Figure [2|c) shows a measurement of T versus the static 
flux detuning /dc and the asymmetry parameter a. The 
static component /dc is used to shift the extrema of the 
waveform relative to the avoided crossing. This is illus- 
trated in Figs. |2|d-g), where we sketch the qubit energy 
bands and the flux excursion for three values of /dc- 

For positive /d c , only the negative dip of the bihar- 
monic signal reaches the avoided crossing. This section of 
the wave will generate two Landau-Zener transitions just 
like in Fig.[TJa), resulting in a transition rate T that oscil- 
lates as a function of the interference phase <p picked up 
between the scattering events [shaded area in Figs.[2|e)]. 
At the maximum of T marked by I in Fig. |2|c), the drive 
just barely reaches the avoided crossing, giving a phase 
pickup of ip = 7r. As /dc is decreased, the flux excur- 
sion between scattering events gets longer and the phase 
pickup increases. The next maximum (point II) corre- 
sponds to (p = 37r, and successive maxima occur when- 
ever the interference phase equals odd integers of tt. In 
this simplified picture, we do not consider the phase ac- 
cumulation during the positive part of the flux sweep, 
which has a smaller influence on the resulting transition 
rate due to the qubit' s comparatively short dephasing 
time (T « 10 ns). 

For negative /dc the positive part of the biharmonic 
signal reaches the avoided crossing, as illustrated in 
Fig. [2|g). Depending on parameters, the waveform may 
drive the qubit through the avoided crossing up to four 
times per cycle, giving three interference phases. The 
phase accumulation and the interference conditions vary 
strongly with waveform shape, and this is the origin of 
the rich fluctuations in T around point III in Fig. fflc). 
Coming back to the trajectories discussed in Figs, yjd- 
e), we can interpret the fluctuations in T as interfer- 
ence from all possible paths generated by the scatter- 
ing events. Note that the drive waveform is periodic 
in time, so the number of possible scattering events is 
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FIG. 3: Fluctuations in the transition rate of a driven qubit. 
(a) Measured transition rate vs flux detuning and drive asym- 
metry a. (b) Transition rate averaged from —4 m$o to m$o- 
The data does not show any dependence on the drive asym- 
metry, (c) Standard deviation of the transition rate. The 
fluctuations have a sharp peak at a = 0. 

not restricted to the four events discussed in Figs, [ljd- 
e). Rather, interference occurs between all trajectories 
that retain phase-coherence, and the maximal number of 
paths is ultimately set by the coherence time of the qubit. 

Figure [3ja) shows a magnification of the region around 
point III in Fig. [2|c). Despite involving only a few scat- 
tering events, the plot shows a rich interference pattern 
as a function of /d c and a. To make the connection to 
charge transport in mesoscopic systems, we identify the 
qubit's transition rate T with the electric conductance, 
whereas the time-symmetry breaking parameter a cor- 
responds to the magnetic field. We calculate the mean 
and the standard deviation of T by averaging over /d c in 
the range -4m$o < /dc < 0m$o. Since the parame- 
ter /dc controls the timing of the scattering events, this 
averaging effectively corresponds to averaging over dif- 
ferent scattering configurations. In mesoscopic systems, 
the averaging over different scattering configurations is 
typically done by tuning a gate voltage or the in-plane 
magnetic field. 

Figures [3jb) and|3|c) show the average transition rate 
(r) and the standard deviation of T, respectively, as a 
function of drive waveform asymmetry. We note that 
(r) is independent of a. This is in contrast to mesoscopic 
systems, which typically show a dip in the conductance 
due to weak localization in the presence of time-reversal 
symmetry. Weak localization is thus not present in our 
experiments. On the other hand, the sharp peak in the 
standard deviation of T at a = [Fig. [3^c)] corresponds 
to the UCF-peak of phase-coherent electron transport. 
The existence of UCF-like features without weak local- 
ization is related to our particular implementation of the 
driving protocol. For a biharmonic drive signal, the av- 
erage qubit transition rate is given by a sum of Bessel 
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FIG. 4: Simulation results, (a) Simulated transition rate F 
for qubit parameters T2 = 10 ns and A/h = 19 MHz. (b) 
Standard deviation of the simulated transition rates for dif- 
ferent values of T2, extracted from simulation results similar 
to those shown in Fig. [4^a). Note that the curves are not 
offset from each other. 

functions, which satisfy a sum rule p~8| [19] that holds 
regardless of waveform asymmetry. We calculate the a- 
independent (r) = 0.15Mrad/s [22], in agreement with 
the results of Fig. |3|c). We expect weak localization to 
be present when driving the qubit with a more arbitrary 
waveform, in which case the average transition rate will 
not be constrained by a sum rule. 

Measurements of weak localization and universal con- 
ductance fluctuations are normally used to determine 
coherence lengths in two-dimensional systems. For our 
driven qubit, where the interfering trajectories are in 
phase space rather than real space, we can employ the 
technique to investigate the coherence time T2 (~ T^) of 
the device. To this end, we have simulated the qubit's 
transition rate numerically [19] for different values of T2 
under the drive defined by Eq. 0. Figure [4^ a) shows 
the result of such a simulation, with qubit parameters 
A/h= 19 MHz (h is Planck's constant), T\ = 20 /is, and 
T2 = 10 ns, in striking agreement with the measured data 
in Fig. [3(a). 

In Fig. [4(b), we plot the standard deviation of the tran- 
sition rate, calculated by repeating the simulation for dif- 
ferent values of T2. The peak around zero asymmetry is 
very broad for small values of T2, but becomes sharper 
and stronger as T 2 is increased. Features due to weak lo- 
calization and UCF in mesoscopic systems show the same 
behavior. For long T2, there are structures appearing in 
Fig. [4(b), away from the main peak at a = 0, which 
are an effect of the small number of possible trajectories 
in our setup. In mesoscopic systems, the large number 
of scattering possibilities average out all configuration- 
specific structures, and leaves only features related to 
the breaking of time-reversal symmetry. We note that 
the simulations give good agreement with the data for 
T2 = 10 — 20 ns, consistent with previous results [19]. 



To conclude, we have investigated effects due to time- 
reversal symmetry in a driven two-level system where 
both the scattering events and the time-reversal sym- 
metry are imposed by the driving protocol. This setup 
allows us to perform experiments in the regime of few 
scatterers, with a finite and controllable number of pos- 
sible trajectories. We find that effects similar to univer- 
sal conductance fluctuations, normally associated with 
chaotic systems with a large number of scattering events, 
persist even in the few-scatterer limit. In addition, our 
work shows an example of how a well-controlled, well- 
understood device like a qubit can be used to simulate 
more complex quantum systems. 

We thank T. Orlando, M. Rudner and L. Levitov for 
helpful discussions. 
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